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A  3D  mathematical  model  of  an  ordered  nanostructured  cathode  catalyst  layer  (CCL)  has  been  developed 
for  proton  exchange  membrane  (PEM)  fuel  cells.  In  an  ordered  nanostructured  CCL,  carbon  nanotubes 
(CNTs)  are  used  as  support  material  for  Pt  catalyst,  upon  which  a  thin  layer  of  proton-conducting  polymer 
(Nation)  is  deposited,  which  are  then  aligned  along  the  main  transport  direction  (perpendicular  to  the 
membrane)  of  various  species.  The  model  considers  all  the  relevant  processes  in  different  phases  of  an 
ordered  nanostructured  CCL.  In  addition,  the  effect  of  Knudsen  diffusion  is  accounted  in  the  model.  The 
model  can  predict  not  only  the  performance  of  an  ordered  nanostructured  CCL  at  various  operating  and 
design  conditions  but  also  can  predict  the  distributions  of  various  fields  in  different  phases  of  an  ordered 
nanostructured  CCL.  The  predicted  nanostructured  CCL  performance  with  estimated  membrane  overpo¬ 
tential  is  validated  with  measured  data  found  in  the  literature,  and  a  good  agreement  is  obtained  between 
the  model  prediction  and  measured  result.  Moreover,  a  parametric  study  is  conducted  to  investigate  the 
effect  of  key  design  parameters  on  the  performance  of  an  ordered  nanostructured  CCL.  In  the  absence 
of  liquid  water,  it  is  found  that  oxygen  diffusion  in  the  pore  phase  is  not  the  limiting  factor  for  the  per¬ 
formance  of  an  ordered  nanostructured  CCL,  owing  to  its  parallel  gas  pores  and  high  porosity.  However, 
the  transport  of  dissolved  oxygen  through  the  Nation  phase  has  a  significant  effect  on  the  performance 
of  an  ordered  nanostructured  CCL.  Further,  it  is  found  that  increasing  the  spacing  between  CNTs  results 
in  a  considerable  drop  in  the  performance  of  an  ordered  nanostructured  CCL  at  the  base  case  conditions 
considered  in  the  simulation. 

©  201 1  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

In  spite  of  significant  improvements  in  the  performance  of 
PEM  fuel  cells  over  the  last  decade,  cost  and  durability  are  still 
the  most  critical  barriers  hindering  their  commercialization  [1,2]. 
Cost  reduction  can  be  realized  by  more  efficient  use  of  materials 
in  the  electrodes  of  PEM  fuel  cells  and  economy-of-scale  in  pro¬ 
duction  [3].  For  instance,  maximizing  the  platinum  (Pt)  catalyst 
utilization  in  the  catalyst  layers  (CLs)  of  PEM  fuel  cells,  results  in 
reducing  the  Pt  loading,  which  in  turn  helps  in  cost  savings.  On 
the  other  hand,  durability  of  PEM  fuel  cells  depends  on  all  of  their 
membrane-electrode-assembly  (MEA)  components.  However,  the 
major  contributors  to  the  excessive  degradation  of  stack  voltage 
are  the  Pt  surface  area  loss  due  to  carbon  corrosion  and  Pt  dissolu¬ 
tion  (aggregation)  in  CLs  [4].  Moreover,  of  the  two  CLs,  CCL  is  the 
source  of  largest  potential  loss  in  PEM  fuel  cells,  owing  to  sluggish 
kinetics  of  oxygen  reduction  reaction  (ORR)  and  the  mass  transport 
limitations  due  to  the  formation  and  transport  of  liquid  water. 
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The  conventional  state-of-the-art  PEM  fuel  cell  CCL  contains  cat¬ 
alyst,  electron  conductor,  proton  conductor  and  pores.  Typically, 
carbon  black  is  used  as  the  electron  conductor  on  which  2-3  nm 
sized  Pt  catalyst  particles  are  finely  dispersed.  It  is  usually  mixed 
with  proton  conducting  polymer  like  Nation  and  a  solvent  to  form 
an  ink.  Pores  are  formed  during  drying  of  an  ink,  after  the  solid 
particles  reached  their  maximum  packing  density  [3].  Despite  its 
extensive  use,  carbon  black  in  the  conventional  CCL  is  reported  to 
undergo  electrochemical  oxidation  to  surface  oxides,  and  finally  to 
C02.  Consequently,  Pt  particle  agglomerates  are  formed  resulting  in 
the  loss  of  electrochemical  surface  area  (ECSA)  [5-1 3  ].  Additionally, 
the  morphology  of  the  conventional  CCL  studied  through  high  reso¬ 
lution  scanning  electron  microscopy  (HR-SEM),  reveals  that  a  large 
portion  of  the  Pt  catalyst  is  located  deep  inside  the  agglomerates, 
and  is  hardly  accessible  to  reactant  oxygen  leading  to  inefficient  Pt 
utilization  [3,14,15].  Further,  there  will  be  increased  mass  transfer 
resistance  for  the  reactant  gas  to  reach  to  the  Pt  particles  located 
deep  inside  the  agglomerates,  resulting  in  high  concentration  over¬ 
potential. 

In  order  to  overcome  the  issues  associated  with  conventional 
CCL,  such  as  carbon  corrosion  leading  to  loss  of  ECSA,  low  Pt  utiliza¬ 
tion  resulting  in  wastage  of  expensive  Pt,  and  transport  limitations 
resulting  in  overall  performance  loss,  considerable  efforts  have 
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Nomenclature 

A; 

catalyst  surface  area  per  unit  mass  of  the  catalyst, 
m2  g-1 

Av 

reactive  surface  area  per  unit  volume,  m2  m-3 

c 

concentration  of  the  mixture  in  the  pore  phase, 
mole  m3 

Co 

oxygen  concentration  at  the  GDL/CCL  interface, 
mole  m3 

CH+ 

fixed  charge  concentration  in  the  membrane, 
mole  m3 

C°2  -g 

gaseous  oxygen  concentration  in  the  pore  phase, 
mole  m3 

c0  2~m 

dissolved  oxygen  concentration,  mole  m-3 

C0  2, ref 

reference  oxygen  concentration,  mole  m-3 

dP 

mean  pore  diameter  of  the  pore  equivalent  to  spac¬ 
ing  between  two  CNTs,  nm 

dcNTs 

distance  or  spacing  between  two  CNTs,  nm 

D,j 

ordinary  diffusion  coefficient,  m2  s-1 

Dl<n,i 

Knudsen  diffusion  coefficient  of  species  i,  m2  s-1 

D02-m 

oxygen  diffusion  coefficient  in  the  Nation  phase, 
cm2  s-1 

(Pij 

combined  diffusion  coefficient,  m2  s-1 

KDeff 

Pij 

Er 

effective  combined  diffusion  coefficient,  m2  s-1 

thermodynamic  reversible  cell  potential,  V 

F 

Faraday’s  constant,  96,487  C  mole-1 

"o2_m 

Flenry’s  constant  for  oxygen  dissolution  in  the 
Nation  film,  J  mole-1 

Jm 

protonic  current  density  in  the  Nation  phase,  Am-2 

href 

reference  exchange  current  density  for  ORR,  Am-2 

k 

constant,  3.16  x  10-8  Pam2  s-1 

I<e 

electrokinetic  permeability  of  the  membrane,  m2 

KP 

hydraulic  permeability  of  the  membrane,  m2 

Ecm 

length  of  CNT  or  CCL  thickness,  p,m 

mPt 

catalyst  loading  per  unit  area  of  the  cathode, 
mg  cm-2 

Mi 

molecular  weight  of  species  i,  kg  mol-1 

Ni 

molar  flux  of  species  i  in  the  pore  phase, 
mole  m  2  s-1 

No2_m 

molar  flux  of  dissolved  oxygen  in  the  Nation  phase, 
mole  m  2  s-1 

P 

pressure,  Pa 

Ph2 

partial  pressure  of  hydrogen,  atm 

Po2 

partial  pressure  of  oxygen,  atm 

rCNT 

radius  of  CNT,  nm 

rb 

interface  between  the  pore  phase  and  the  Nation 
phase 

R 

universal  gas  constant,  8.3143  J  mole-1 1<-1 

tm 

membrane  thickness,  |jim 

^Nafion 

nation  film  thickness,  nm 

T 

temperature,  K 

X 

x-direction,  m 

mole  fraction  of  species  i  in  the  pore  phase 

y 

y-direction,  m 

Wf 

mass  fraction  of  species  i  in  the  pore  phase 

z 

z-direction,  m 

Greek 

Ota 

anode  charge  transfer  coefficient 

Otc 

cathode  charge  transfer  coefficient 

Apa-c 

pressure  difference  between  the  anode  and  the 
cathode,  Pa 

y 

reaction  order  for  ORR 

£g 

porosity  or  volume  fraction  of  pore  phase 

£m  volume  fraction  of  Nation  phase 

Tjact  cathode  activation  overpotential,  V 

rjohm  cathode  ohmic  overpotential,  V 

tjq  applied  cathode  overpotential  at  the  membrane 

interface,  V 

ijm  membrane  ohmic  overpotential,  V 

k  intrinsic  proton  conductivity  of  Nation,  S  nrr 1 

Keff  effective  proton  conductivity  in  the  Nation  phase, 
Sm-1 

/xh2o  viscosity  of  the  liquid  water,  kg  nrr1  s-1 

Vi  molar  diffusion  volume  of  species  i,  m3  mole-1 

0m  protonic  phase  potential  in  the  Nation  phase,  V 
0s  electronic  phase  potential  in  the  electron¬ 
conducting  phase  (CNT),  V 


been  dedicated  to  develop  novel  catalyst  support  material  for  CCL, 
and  the  resulting  CCL  is  often  referred  to  as  ordered  nanostructured 
or  oriented  CCL  structure  [3,4,9,16-32].  In  an  ordered  nanostruc¬ 
tured  CCL,  carbon  nanotubes  (CNTs)  are  being  used  as  support 
material  for  Pt  catalyst  owing  to  their  unique  features,  which 
include  superior  mechanical  and  electrical  properties,  high  aspect 
ratio  and  enhanced  mass  transport  capability  [16,18,19,33-38].  On 
the  surface  of  these  CNTs,  also  referred  to  as  electron  conductors,  Pt 
particles  of  approximately  2  nm  are  finely  dispersed.  A  thin  layer  of 
Nation  is  then  deposited  on  the  Pt-coated  CNTs  and  aligned  along 
the  main  transport  direction  (perpendicular  to  the  membrane)  of 
electron,  proton  and  reactant  gas  to  form  an  ordered  or  oriented 
nanostructured  CCL. 

Generally,  it  is  believed  that  an  oriented  nanostructured  CCL 
offers  much  improved  performance  than  a  non-oriented  or  disor¬ 
dered  nanostructured  CCL  structure  because  of  its  following  salient 
features.  First,  the  inherent  advantage  of  an  oriented  nanostruc¬ 
tured  CCL,  resulting  in  straight  pores,  would  enhance  the  transport 
of  oxygen.  Second,  the  Pt  particles  are  located  at  the  interface  of 
the  electron  conductor  (CNTs)  and  the  proton  conductor  (Nation) 
would  result  in  better  Pt  utilization  in  an  oriented  nanostructured 
CCL.  Third,  the  electronic  conductivity  of  CNTs  is  much  higher  along 
the  tubes  than  across  the  tubes;  as  a  result,  there  is  negligible  loss  in 
potential  when  electrons  transfer  along  the  tubes  [39,40].  Fourth, 
an  oriented  nanostructured  CCL  may  exhibit  superhydrophobicity, 
which  can  facilitate  greater  water  removal  within  the  CCL,  thereby 
improving  mass  transport  [18]. 

Mathematical  modeling  is  an  essential  tool  in  understanding 
the  processes  occurring  in  different  phases  of  an  oriented  nanos¬ 
tructured  CCL  and  helps  in  gaining  insights  about  the  processes 
influencing  the  CCL  performance.  In  the  recent  past,  there  were  few 
attempts  to  model  an  ordered  nanostructured  CCL.  Du  et  al.  [41 ,42] 
were  few  of  the  earliest  researchers  to  model  an  ordered  nanostruc¬ 
tured  CCL  having  straight  pores.  Later,  Chisaka  and  Daiguji  [43]  also 
modeled  an  ordered  nanostructured  CCL  having  triangular  orienta¬ 
tion  of  CNTs.  These  models  [41  -43  ]  were  1 D  models  and  considered 
mostly  the  radial  diffusion  of  oxygen  through  the  Nation  film.  More¬ 
over,  Chisaka  and  Daiguji  [43]  did  not  consider  the  axial  gradient 
of  oxygen  concentration  induced  along  the  pores.  More  recently, 
Rao  and  Xing  [44]  developed  a  2D  model  for  an  ordered  nanos¬ 
tructured  CCL  having  square  orientation  of  CNTs.  They  assumed  a 
thick  Nation  film  (1  p,m)  on  the  surface  of  CNTs,  which  is  inconsis¬ 
tent  with  the  reported  values  for  an  ordered  nanostructured  CCL 
and  also  neglected  the  proton  transport  in  the  Nation  film.  Further, 
none  of  the  models  [41  -44]  were  validated  with  the  measured  data 
available  in  the  literature. 

The  objective  of  the  present  study  is  to  develop  a  3D  model  of 
an  ordered  nanostructured  CCL  for  PEM  fuel  cells,  which  considers 
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Membrane 


Fig.  1.  Schematic  representation  of  an  ordered  nanostructured  CCL. 


all  the  relevant  processes  in  different  phases  of  the  nanostructured 
CCL  and  simulates  its  performance  using  realistic  design  conditions. 
Moreover,  the  other  objectives  are  to  validate  the  predicted  nanos¬ 
tructured  CCL  performance  with  measured  data  available  in  the 
literature  and  to  assess  the  effects  of  key  morphological  parameters 
on  the  performance  of  an  ordered  nanostructured  CCL. 

2.  Model  formulation 

A  schematic  representation  of  an  ordered  nanostructured  CCL 
is  shown  in  Fig.  1.  It  consists  of  oriented  CNTs,  upon  which  Pt  par¬ 
ticles  are  uniformly  distributed,  and  the  catalyzed  CNTs  are  finally 
covered  with  thin  films  of  Nation.  Oxygen  in  the  multi-component 
mixture  diffuses  through  the  straight  pores  in  between  the  CNTs, 
and  then  diffuses  radially  through  the  Nation  film  to  reach  the  reac¬ 
tion  sites  on  the  surface  of  CNTs. 

In  an  ordered  nanostructured  CCL,  CNTs  can  be  oriented  in  any 
particular  fashion.  However,  in  the  present  study,  a  square  grid  ori¬ 
entation  of  CNTs  is  considered,  which  allows  using  symmetry  in  the 
simulation.  Fig.  2  illustrates  the  cross-sectional  view  of  an  ordered 
nanostructured  CCL;  whereas,  Fig.  3  shows  the  3D  computational 
domain.  The  dashed  line  square  in  Fig.  2  shows  the  symmetric  cross- 
sectional  portion  of  an  ordered  nanostructured  CCL  in  they-z  plane. 
The  computational  domain  includes  the  quarter  portion  of  CNT,  the 
Nation  thickness,  and  half-spacing  between  two  CNTs  representing 
the  pore-phase.  Moreover,  the  x-direction  represents  the  length  of 
CNT  or  CCL  thickness.  Oxygen  in  the  multi-component  mixture  dif¬ 
fuses  into  the  CCL  through  the  gas  diffusion  layer  (GDL)  on  the 
left  boundary  (i.e.  x  =  0)  and  protons  migrate  into  the  CCL  from 
the  membrane  side  on  the  right  boundary  (i.e.  x  =  LCnt).  Oxygen 
then  dissolves  in  the  Nation  film  and  is  consumed  at  the  reaction 
boundary  (rCNT)  along  with  protons  and  electrons  due  to  the  oxygen 
reduction  reaction  (ORR). 

The  cell  is  assumed  to  operate  under  steady  state  condition.  The 
temperature  and  pressure  are  assumed  to  be  uniform  in  the  CCL. 
The  convective  flux  is  assumed  to  be  negligible  when  compared 
to  the  diffusive  flux  of  species.  The  multi-component  mixture  is 
approximated  as  ideal  gas  with  negligible  viscous,  Soret,  Dufour 
and  gravity  effects.  The  product  water  is  assumed  to  be  in  the  gas 
phase  and  diffuses  out  of  the  CCL  sufficiently  fast,  implying  there  is 
no  liquid  water  present  in  the  pore  phase  of  an  ordered  nanostruc¬ 
tured  CCL.  The  proton  conductivity  in  the  Nation  phase  is  assumed 
to  be  constant,  and  due  to  sufficiently  high  electronic  conductivity 
of  CNTs,  the  electronic  phase  potential  is  assumed  to  be  negligi¬ 
ble  compared  to  the  protonic  phase  potential  in  the  CCL.  Finally, 
the  Pt  catalyst  is  assumed  to  be  uniformly  distributed  over  the  CNT 
surface,  implying  ORR  occurs  at  the  interface  between  electron¬ 


Nafion  film  thickness 

(tNafion)  Spacing  (dCNis) 


conducting  phase  (CNT)  and  proton-conducting  phase  (Nation), 
also  referred  to  as  the  reaction  boundary. 

Based  on  the  above  stated  assumptions,  the  governing  equa¬ 
tions  for  the  transport  of  multi-component  mixture  in  the  pore 
phase,  dissolved  oxygen  transport  in  the  Nation  phase,  and  proton 
transport  in  the  Nation  phase  are  described  as  follows. 

Multi-component  mixture  transport :  V  •  7?/  =  0  (1) 

Dissolved  oxygen  transport :  V  •  ~^o2_m  =  0  (2) 


x  =  o 

(GDL  Interface) 


Fig.  3.  3D  computational  domain  of  an  ordered  nanostructured  CCL 
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Proton  transport :  V  Jm  =  0  (3) 

where  Nt  is  the  species  diffusion  flux  in  the  pore  phase,  No2_m  is  the 
diffusion  flux  of  dissolved  oxygen  in  the  Nation  phase,  and  Jm  is  the 
protonic  current  density  in  the  Nation  phase. 

The  diffusion  flux  in  the  pore  phase  can  be  determined  using  the 
modified  Stefan-Maxwell  equations  incorporating  Knudsen  diffu¬ 
sion  for  multi-component  systems  involving  n  species,  expressed 
as  [45]: 


VXj  = 


—  xj~^  i)  + 


A 

DI<n,i 


(4) 


where  c  is  the  concentration  of  the  mixture  in  the  pore  phase,  Dy 
is  the  ordinary  diffusion  coefficient  of  species  i  in  j,  x,-  is  the  mole 
fraction  of  species  i  in  the  pore  phase,  and  7? f  is  the  diffusion  flux  of 
species  i  in  the  pore  phase,  and  DKn  i  is  the  Knudsen  diffusion  coeffi¬ 
cient  of  species  i,  which  is  calculated  using  the  following  expression 
[46-48]: 


Di(n,i  — 


(5) 


CNTs  are  uniformly  distributed  and  form  a  straight  pore  network 
along  the  axial  direction,  the  effective  properties  are  corrected  only 
with  volume  fractions  of  pore  or  Nation. 

Expressing  the  dissolved  oxygen  flux  and  protonic  current  den¬ 
sity  in  the  Nation  phase  in  terms  of  dissolved  oxygen  concentration 
and  protonic  phase  potential,  respectively,  Eqs.  (2)  and  (3)  become 

Dissolved  oxygen  transport :  V  •  (Dq2  m  c02_m )  =  0  (11) 

Proton  transport :  V  •  (/ceff0m)  =  0  (12) 

where  Dq2  m  is  oxygen  diffusion  coefficient  in  the  Nation  phase, 
Co2_m  and  0m  are  the  dissolved  oxygen  concentration  and  protonic 
phase  potential  in  the  Nation  phase,  respectively,  and  /ceff  is  the 
effective  proton  conductivity  in  the  Nation  phase  of  an  ordered 
nanostructured  CCL,  defined  as 

/ceff  =  emK  (13) 

where  sm  is  the  volume  fraction  of  Nation  in  an  ordered  nanostruc¬ 
tured  CCL  and  k  is  the  intrinsic  proton  conductivity  of  Nation. 

The  oxygen  diffusion  coefficient  in  the  Nation  phase  (Do2_m )  is 
estimated  as  [50]: 


where  dp  is  the  mean  pore  diameter  equivalent  to  mean  distance 
or  spacing  between  two  CNTs  (dcNTs)  in  an  ordered  nanostruc¬ 
tured  CCL,  R  is  the  universal  gas  constant  in  J  mole-1  K1,  T  is  the 
temperature  in  K,  and  M*  is  the  molecular  weight  of  species  i  in 
kg  mole-1 . 

The  binary  or  ordinary  diffusion  coefficients  (Dy)  can  be  calcu¬ 
lated  using  an  empirical  relation  based  on  the  kinetic  theory  of 
gases,  expressed  as  [49]: 


Du  =  k- 


jl.75 


p(v,/3  +  v,/3 


1  1 

Mi  +  Mi 


1 

2 


(6) 


where  k  is  a  constant  with  the  value  3.16x1 0-8  Pa  m2  s-1 ,  T  is  the 
temperature  in  K,  p  denotes  the  pressure  in  Pa,  V;  is  the  molar  dif¬ 
fusion  volume  of  species  i  in  m3  mole-1,  and  M,-  is  the  molecular 
weight  of  species  i  in  kg  mole-1. 

Mathematically,  the  two  terms  on  the  right  hand  side  of  Eq.  (4) 
can  be  combined  to  obtain  the  modified  Stefan-Maxwell  model  for 
multi-component  diffusion  in  the  pore  phase  of  an  ordered  nanos¬ 
tructured  CCL,  and  can  be  expressed  as 


VXj  = 


-xfli) 


(7) 


where  Py  is  the  combined  diffusion  coefficient,  defined  as 


DO2_m=3.1xl0-3exp(-^)  (14) 

where  T  is  in  K,  and  Dnn  is  in  cm2  s-1 . 

The  cathode  potential  in  the  model  is  defined  as 


Cathode  potential  =  Er  -  r\act  -  dohm 


(15) 


where  Er  is  the  thermodynamic  reversible  cell  potential,  pact  is  the 
potential  loss  due  to  resistance  to  electrochemical  reaction  and 
mass  transfer  limitations  in  the  CCL,  also  referred  here  as  cathode 
activation  overpotential,  and  r/ohm  is  the  potential  loss  due  to  proton 
transport  in  the  CCL,  also  referred  to  as  cathode  ohmic  overpoten¬ 
tial,  which  can  be  found  using  the  Ohm’s  law.  The  thermodynamic 
reversible  cell  potential  (Er)  as  a  function  of  operating  temperature, 
pressure  and  reactant  concentrations,  can  be  expressed  as  [50]: 

Er  =  1.23  -  0.9  X  10_3(T  -  298)  +  —  ln(p32Po2)  (16) 

where  pu2  and  po2  are  the  partial  pressures  (in  atm)  of  hydrogen 
and  oxygen,  respectively. 

The  dissolved  oxygen  concentration  at  the  interface  between  the 
pore  phase  and  the  Nation  phase  is  related  to  the  gaseous  oxygen 
concentration  using  the  Henry’s  law,  expressed  as  [50,51  ]: 


C02-mlrh  - 


u2 -mlrfa  LJ  C02-g\rb 

1102 -m 


(17) 
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\Dij  +  DKn  j  J 
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In  terms  of  effective  diffusion  coefficient,  modified 
Stefan-Maxwell  model  is  expressed  as 


n 


where  p?.ff  is  the  combined  effective  diffusion  coefficient  in  the  pore 
phase,  defined  as 

pf  =  egpij  (10) 

where  sg  is  pore  volume  fraction  or  porosity  of  an  ordered  nanos¬ 
tructured  CCL.  It  is  worthwhile  to  note  that  the  effective  properties 
are  not  corrected  for  porosity  and  tortuosity  using  the  Bruggeman 
correlation.  However,  in  an  ordered  nanostructured  CCL,  since  the 


where  c02 _g  |r&  is  the  gaseous  oxygen  concentration  at  the  interface 
of  the  Nation  phase  and  the  pore  phase  and  Hq2  m  is  the  Henry’s  con¬ 
stant  for  oxygen  dissolution  in  the  Nation  film,  which  is  calculated 
using  the  empirical  correlation  [52,53]: 

Ho2_m=  0.1013  exp  (l4.1-^)  (18) 

where  the  temperature  T  is  in  K,  and  Hq2 _m  is  in  Pam3  mole-1  or 
Jmole-1. 

Similarly,  at  the  interface  between  the  Nation  phase  and  CNT 
phase  of  an  ordered  nanostructured  CCL,  ORR  takes  place  and  its 
rate  expression  is  given  by  the  Butler-Volmer  equation,  expressed 
as 


Rc  —  AJo,re/ 


C°2-m  lrCNT 
C°2  ,ref 


(19) 
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Table  1 

Governing  equations. 


Table  2 

Boundary  conditions. 


Multi-component  mixture  transport 

V • T$i  =  0 

n 

Locations 

-cVx,  = 

‘  ^  5  ij 

J=l,i 

x  =  0 

Dissolved  oxygen  transport 

V  ■  (Do2  mCo2  m)  =  0 

Proton  transport 

V.(Keff0m)  =  0 

where  Av  is  the  reactive  surface  area  density,  J0, re/  is  the  reference 
exchange  current  density  measured  at  a  reference  oxygen  concen¬ 
tration  (co2re/),  co2_mlrCNT  the  dissolved  oxygen  concentration 
at  the  reaction  interface,  y  is  the  reaction  order,  which  is  usually 
equal  to  one,  ac  and  aa  are  the  cathodic  and  anodic  transfer  coeffi¬ 
cients,  respectively,  and  rjact  is  the  cathode  activation  overpotential, 
defined  as 

L] act  =  4>s  ~  0m  (20) 
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y  and  z  external  boundaries 


Tcnt 


rb 


Boundary  conditions 


w,-  —  specified 


9c( 


O  2-m 


dx 


=  0 


d(pm 

Jx 

dWf 

§X 

dc( 


=  0 


=  0 


o  2-m 


=  0 


dx 

rj  o  =  specified 
7$i  ■  It  =  0 


o2- 


m  -  It  =  °  > 

It  =  0  J 

RcLcm 


No  2_m  =  - 
Jm  =  -Pc  Lem 

CO; 

c°2-  =  Ho. 
Jm  =  0 


4  F 


2  -g 


where  0S  is  the  electronic  phase  potential  in  the  electron¬ 
conducting  phase  (CNT)  of  an  ordered  nanostructured  CCL. 

The  reactive  surface  area  density  is  given  by  [54]: 


A 


V  — 


mptAs 

Lem 


(21) 


where  mpt  is  the  catalyst  loading  per  unit  area  of  the  cathode,  As  is 
the  catalyst  surface  area  per  unit  mass  of  the  catalyst,  and  LCnt  is 
the  length  of  the  CNT  or  thickness  of  the  CCL  in  an  ordered  nanos¬ 
tructured  CCL 

The  reference  exchange  current  density  (Jo, ref)  f°r  ORR  is 
obtained  using  an  empirical  correlation  based  on  the  experimental 
data  of  Parthasarathy  et  al.  [55],  expressed  as 


logC/o,re/)  =  3.507 


4001 

T 


(22) 


where  J0,re/  is  in  A  cm-2  and  T  is  in  K. 

Finally,  the  governing  equations  in  an  ordered  nanostructured 
CCL  are  summarized  in  Table  1. 


3.  Boundary  conditions 

The  boundary  conditions  required  to  complete  the  model 
formulation  can  be  classified  into  external  and  internal  bound¬ 
ary  conditions.  External  boundary  conditions  are  applied  at  the 
external  boundaries  of  the  computational  domain  shown  in 
Fig.  3.  Whereas,  internal  boundary  conditions  are  applied  at 
the  interfaces  between  the  pore  phase  and  the  Nation  phase, 


and  between  the  CNT  phase  and  the  Nation  phase  of  the  com¬ 
putational  domain,  also  referred  to  as  the  interfacial  boundary 
conditions. 

Since  the  computational  domain  considered  is  a  symmetric 
portion  of  an  ordered  nanostructured  CCL,  symmetric  boundary 
conditions  are  applicable  at  the  external  boundaries  of  y  and  z 
axes.  In  contrary,  at  the  external  boundary  interfacing  with  GDL 
(indicated  by  x  =  0  in  Fig.  3),  mass  fractions  of  multi-component 
mixture  species  are  specified,  protonic  current  density  and  dis¬ 
solved  oxygen  flux  are  set  as  zero.  Similarly,  at  the  external 
boundary  interfacing  with  the  membrane  (indicated  by  x  =  LCnt  in 
Fig.  3),  cathode  overpotential  is  specified,  diffusional  flux  for  multi- 
component  mixture  species  and  dissolved  oxygen  flux  are  set  as 
zero. 

The  locations  at  which  internal  boundary  conditions  are 
required  in  the  model  formulation  are  illustrated  as  rCNT  and  rb. 
At  the  interface  between  the  CNT  phase  and  Nation  phase  (rCNi), 
ORR  takes  place,  also  referred  to  as  the  reaction  boundary,  reac¬ 
tion  boundary  conditions  are  applied  for  the  dissolved  oxygen  and 
proton  transport.  On  the  other  hand,  at  the  interface  between  the 
pore  phase  and  Nation  phase  (rb),  Flenry’s  law  is  applied  to  relate 
the  dissolved  oxygen  concentration  with  the  gaseous  oxygen  con¬ 
centration  in  the  pore  phase,  and  protonic  current  density  is  set 
as  zero.  The  mathematical  form  of  boundary  conditions  applied  is 
tabulated  in  Table  2. 


1“CNT  (nm) 

^CNT  (nm) 


Fig.  4.  Distributions  of  normalized  gaseous  oxygen  concentration  (a)  and  ORR  rate  (b)  along  the  length  of  CNT  in  an  ordered  nanostructured  CCL  at  different  specified 
overpotentials. 
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Table  3 

Parameters  used  for  model  validation. 


Operating  temperature,  T 
Pressure,  p 

Oxygen  concentration  at  the  GDL/CCL 
interface,  Co 

Reference  oxygen  concentration,  co2  ref 
Radius  of  CNT,  rCm 
Thickness  of  the  Nation  film,  tNation 
Distance  between  two  CNTs,  dcNTs 
Length  of  CNT,  LCnt 

Intrinsic  proton  conductivity  of  Nation,  k 
Pt  loading,  mPt 

Specific  surface  area  of  Pt  particles,  As 
Membrane  thickness,  tm 
Fixed  charge  concentration  in  the  membrane, 
cH+ 

Hydraulic  permeability  of  the  membrane,  I<P 
Electrokinetic  permeability  of  the  membrane, 
Ke 


Source:  Li  et  al.  [18]. 


343.0  K 
2.0  atm 

70.13  mole  m-3  (estimated) 

1.0  molem-3 
25  nm 
5  nm 

50  nm  (assumed) 

5  pan 

7.0  5  m-1  [57] 

0.2  mg  cm-2 
100  m2  g_1 
164  p.m  [56] 

1200 molem-3  [56] 

1.8  x  10-18  m2  [56] 

7.18  xlO-20  m2  [56] 


4.  Numerical  implementation 

The  governing  equations  along  with  boundary  conditions  pre¬ 
sented  in  the  above  section  are  solved  using  a  commercial  package, 
commonly  known  as  COMSOL  Multiphysics™.  The  package  uses 
finite  element  method  to  solve  the  governing  partial  differen¬ 
tial  equations  (PDEs)  over  a  defined  geometry.  The  computational 
domain  is  discretized  into  tetrahedral  mesh  and  Lagrangian- 
quadratic  elements  of  second  order  are  used.  Since  the  governing 
PDEs  form  a  set  of  non-linear  equations,  a  non-linear  solver  is  used 
to  obtain  the  solution.  The  solution  is  considered  to  be  a  converged 
when  the  relative  tolerance  reaches  below  10-6. 

Additionally,  a  grid  sensitivity  test  has  been  conducted  to  ensure 
grid  independent  solution.  It  also  ensures  the  adequacy  of  mesh 
elements  used  in  the  3D  computational  domain.  Different  level  of 
mesh  sizes  are  considered  ranging  approximately  from  35,000  to 
0.4  million  degrees  of  freedom  (DOF).  The  grid  sensitivity  results 
are  shown  in  Fig.  4,  which  illustrate  normalized  gaseous  oxygen 
concentration  and  ORR  rate  along  the  length  of  CNT  at  differ¬ 
ent  specified  overpotentials.  It  can  be  observed  from  the  figure 
that  all  the  grid  sizes  exhibit  almost  identical  results,  indicat¬ 
ing  grid  independent  solution.  However,  for  better  accuracy,  the 
results  presented  in  the  subsequent  figures  are  for  higher  grid  level 
(DOF  =  71 104). 

5.  Model  validation 

Model  validation  is  as  important  as  model  development.  It  helps 
in  determining  the  range  of  validity  and  accuracy  of  the  model. 
However,  none  of  the  previous  models  [41  -44]  of  an  ordered  nanos- 
tructured  CCL  attempted  to  validate  or  compare  their  predictions 
with  the  measured  data  available  in  the  literature.  The  parameters 
used  in  the  validation  of  the  present  model  are  listed  in  Table  3. 

Most  of  the  parameters  used  in  the  model  validation  are 
obtained  from  Li  et  al.  [18],  who  reported  the  measured  cell 
performance  data  having  an  ordered  nanostructured  CCL.  Since 
the  present  model  predicts  an  ordered  nanostructured  CCL  per¬ 
formance,  to  compare  with  the  measured  cell  performance,  the 
membrane  ohmic  overpotential  is  accounted  using  an  expression 
presented  by  Baschuk  and  Li  [56]: 

C/m/ k)  +  (FKpCH+  Apa-c I K \  (23) 

1  +  (F2/C£c2+//c/xh2o)  / 

where  tm  is  the  membrane  thickness,  I<p  is  the  hydraulic  perme¬ 
ability  of  the  membrane,  cH+  is  the  fixed  charge  concentration  in 
the  membrane,  A pa-c  is  the  pressure  difference  between  the  anode 


Fig.  5.  Comparison  between  the  predicted  and  measured  cell  performance  having 
an  ordered  nanostructured  CCL. 

and  the  cathode  /xH2o  is  the  viscosity  of  liquid  water,  and  I<E  is  the 
electrokinetic  permeability  of  the  membrane. 

Fig.  5  shows  the  comparison  between  the  predicted  and  mea¬ 
sured  cell  performance  having  an  ordered  nanostructured  CCL.  The 
lines  represent  the  model  predictions,  while  the  symbol  repre¬ 
sents  the  measured  cell  performance  data  obtained  from  Li  et  al. 
[18].  The  solid  line  represents  the  predicted  ordered  nanostruc¬ 
tured  CCL  performance  at  343  K  and  2  atm.  When  the  contribution 
of  membrane  overpotential  is  taking  into  account  in  addition  to 
the  cathode  overpotential  towards  the  total  cell  potential  loss,  the 
predicted  performance  curve  is  shown  by  the  dashed  line,  referred 
to  as  the  ‘anode-free  cell  potential’  in  Fig.  5.  It  can  be  seen  that 
the  predicted  anode-free  cell  performance  exhibits  similar  trend  as 
the  measured  cell  performance  with  a  constant  difference  through¬ 
out  the  current  density  range,  indicating  the  difference  originating 
from  the  potential  difference  between  the  calculated  thermody¬ 
namic  reversible  cell  potential  and  the  measured  open-circuit  cell 
potential  at  zero  load,  also  referred  to  as  the  mixed  electrode  poten¬ 
tial  in  the  literature  [58-60].  It  is  reported  to  be  in  the  order  of 
0.1 -0.2  V,  and  arises  due  to  unavoidable  parasitic  reactions  and 
fuel  crossover  through  the  membrane  [53,60].  When  the  anode- 
free  cell  potential  is  corrected  with  the  mixed  electrode  potential, 
the  resulting  cell  performance  curve  is  shown  by  the  dashdot  line 
in  Fig.  5,  and  referred  to  as  the  anode-free  corrected  cell  potential. 
Further,  it  can  be  seen  that  anode-free  corrected  cell  potential  is  in 
good  agreement  with  the  measured  data. 

6.  Results  and  discussion 

The  developed  model  not  only  predicts  the  distributions  of  var¬ 
ious  fields  in  different  phases  of  an  ordered  nanostructured  CCL 
but  also  predicts  the  ordered  nanostructured  CCL  performance  at 
various  operating  and  design  conditions.  The  base  case  parameters 
used  in  the  simulation  are  listed  in  Table  4. 

Fig.  6  shows  the  distributions  of  gaseous  oxygen  concentration 
in  the  pore  phase  of  an  ordered  nanostructured  CCL  at  various 
applied  overpotentials  at  the  membrane  interface.  The  horizontal- 
axis  represents  the  length  of  CNT,  can  also  be  referred  as  the 
thickness  of  the  nanostructured  CCL,  while  the  vertical-axis  repre¬ 
sents  the  gaseous  oxygen  concentration  (c02_g)  in  the  pore  phase. 
The  effect  of  Knudsen  diffusion  is  not  considered  in  the  base  case 
simulation.  A  more  detailed  examination  of  the  effect  Knudsen 
diffusion  with  CNT  spacing  (dcNTs)  is  presented  in  the  following 


Table  4 

Base  case  parameters  used  in  the  simulation. 
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Operating  temperature,  T 

353.0  K 

Pressure,  p 

1.0  atm 

Molar  diffusion  volume  of  oxygen,  v02 

16.6xl0-6m3  mole-1  [61] 

Molar  diffusion  volume  of  water  vapor,  Vh2o 

12.7  x  10-6  m3  mole-1  [61  ] 

Molar  diffusion  volume  of  nitrogen,  vN2 

17.9  x  10-6  m3  mole-1  [61] 

Mass  fraction  of  oxygen  at  the  GDL/CCL  interface, 

0.1447 

Wq2 

Mass  fraction  of  water  vapor  at  the  GDL/CCL 

0.3789 

interface,  wH2o 

Reference  oxygen  concentration,  Co2  ref 

2.637  molem-3 

Cathode  transfer  coefficient,  ac 

0.58 

Anode  transfer  coefficient,  aa 

1 

Radius  of  CNT,  rCm 

25  nm 

Thickness  of  the  Nation  film,  tNation 

lOnm 

Distance  between  two  CNTs,  dcNTs 

50  nm 

Length  of  CNT,  LCnt 

5  pjn 

Pore  volume  fraction,  eg 

0.739 

Nation  volume  fraction,  em 

0.130 

Intrinsic  proton  conductivity  of  Nation,  k 

7.0  Sm-1 

Pt  loading,  mPt 

0.2  mg  cm-2 

Specific  surface  area  of  Pt  particles,  As 

100  m2  g-1 

section.  It  can  be  seen  that  there  is  negligible  variation  in  the  con¬ 
centration  of  gaseous  oxygen  {co2g ),  as  we  move  along  the  length  of 
CNT  for  T]0  =  0.3  V.  However,  with  the  increase  of  applied  overpoten¬ 
tial  at  the  membrane  interface,  the  gaseous  oxygen  concentration 
(co2_g )  decreases  along  the  length  of  CNT  with  a  substantial  drop  at 
?7o  =  0.7V.  Increasing  the  applied  overpotential  (r]0)  increases  the 
ORR  rate  at  the  reaction  boundary  along  the  length  of  CNT  (can 
be  evident  from  Fig.  7),  resulting  in  the  increased  consumption 
of  dissolved  oxygen,  which  in  turn  decreases  the  concentration  of 
gaseous  oxygen  (c02_  )  in  the  pore  phase.  Moreover,  this  is  due  to 
multi-component  diffusion  in  the  pore  phase,  wherein  the  concen¬ 
tration  of  gaseous  oxygen  depends  on  the  concentration  of  other 
n  - 1  species.  Similarly,  the  distributions  of  dissolved  oxygen  con¬ 
centration  at  the  reaction  boundary  along  the  length  of  CNT  at 
various  applied  overpotentials  is  shown  in  Fig.  8.  The  horizontal- 
axis  represents  the  length  of  CNT,  while  the  vertical-axis  represents 
the  dissolved  oxygen  concentration  (co2_m).  Again,  similar  trends 
can  be  observed  for  the  dissolved  oxygen  concentration  (c02_m) 
with  increasing  applied  overpotential  {rj 0)  at  the  membrane  inter¬ 
face. 

On  the  other  hand,  the  distributions  of  dissolved  oxygen  concen¬ 
tration  (co2_m )  within  the  Nation  film  of  an  ordered  nanostructured 


Fig.  6.  Distributions  of  gaseous  oxygen  concentration  along  the  length  of  CNT  in  an 
ordered  nanostructured  CCL  at  different  specified  overpotentials. 


Fig.  7.  Distributions  of  ORR  rate  at  the  reaction  boundary  along  the  length  of  CNT 
in  an  ordered  nanostructured  CCL  at  different  specified  overpotentials. 

CCL  at  different  applied  overpotentials  is  shown  in  Fig.  9.  The 
horizontal-axis  represents  the  thickness  of  the  Nation  film,  while 
the  vertical-axis  represents  the  dissolved  oxygen  concentration 
(co2  m  )•  It  can  be  seen  that  the  dissolved  oxygen  concentration 
(co2_J  exhibits  linear  variation  along  the  thickness  of  the  Nation 
film  at  all  the  applied  overpotentials  (r]0).  Moreover,  it  can  be  seen 
that  the  dissolved  oxygen  concentration  decreases  as  we  move 
from  right  to  the  left,  indicating  the  consumption  of  dissolved  oxy¬ 
gen  at  the  reaction  boundary  (corresponds  to  tNafion  =  0  in  Fig.  9). 
Further,  it  is  worth  mentioning  that  the  concentration  of  dissolved 
oxygen  at  the  right  (corresponds  to  tNafion  =  10  nm  in  Fig.  9)  is  gov¬ 
erned  by  the  Henry’s  law,  indicating  the  interface  between  the 
Nation  phase  and  the  pore  phase. 

7.  Effect  of  Knudsen  diffusion 

Knudsen  diffusion  becomes  important  when  mean  pore  diam¬ 
eter  in  porous  media  (in  other  words,  mean  distance  or  spacing 
between  two  CNTs  in  an  ordered  nanostructured  CCL)  is  compara¬ 
ble  to  the  mean  free  path  of  the  gas  molecule.  According  to  Jennings 


Fig.  8.  Distributions  of  dissolved  oxygen  concentration  at  the  reaction  boundary 
along  the  length  of  CNT  in  an  ordered  nanostructured  CCL  at  different  specified 
overpotentials. 
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Fig.  9.  Distributions  of  dissolved  oxygen  concentration  within  the  Nafion  film  of  an 
ordered  nanostructured  CCL  at  different  specified  overpotentials. 

[62]  expression,  the  mean  free  path  of  air  molecule  at  353  K  and 
1  atm  is  about  80  nm,  which  is  comparable  to  the  mean  distance 
or  spacing  between  two  CNTs  (dCNTs)  at  the  base  case  simulation 
conditions,  indicating  the  diffusion  to  be  in  the  transition  regime. 
Fig.  10  shows  the  effect  of  Knudsen  diffusion  on  the  distribution 
of  gaseous  oxygen  concentration  in  the  pore  phase  of  an  ordered 
nanostructured  CCL  at  an  applied  overpotential  of  0.7  V  using  the 
base  case  parameters  listed  in  Table  4.  Again,  the  horizontal-axis 
represents  the  length  of  CNT,  while  the  vertical-axis  represents 
the  gaseous  oxygen  concentration  (c02_g)  in  the  pore  phase.  Fur¬ 
ther,  the  solid  line  represents  the  concentration  of  gaseous  oxygen 
when  Knudsen  diffusion  is  not  taken  into  account  with  the  ordinary 
diffusion;  whereas,  the  dashed  line  represents  the  concentration 
of  gaseous  oxygen  (c02_g)  when  Knudsen  diffusion  is  taken  into 
account  with  the  ordinary  diffusion  in  the  pore  phase.  It  can  be 
seen  that  the  difference  in  the  concentration  of  oxygen  between 
the  Knudsen  and  non-Knudsen  diffusion  increases  as  we  move 
along  the  length  of  the  CNT.  However,  the  average  drop  in  oxygen 
concentration  along  the  length  of  CNT  between  the  Knudsen  and 


Fig.  10.  Effect  of  Knudsen  diffusion  on  the  gaseous  oxygen  concentration  in  the  pore 
phase  of  an  ordered  nanostructured  CCL  at  r]0  =  0.7  V  using  the  base  case  parameters 
listed  in  Table  4. 


Fig.  11.  Effect  of  Knudsen  diffusion  on  the  performance  of  an  ordered  nanostruc¬ 
tured  CCL  at  the  base  case  conditions  listed  in  Table  4. 

non-Knudsen  diffusion  is  of  the  same  order  of  magnitude,  resulting 
in  a  negligible  effect  on  the  performance  of  an  ordered  nanostruc¬ 
tured  CCL  at  the  base  case  conditions,  which  can  be  evident  from 
Fig.  11.  Additionally,  the  Knudsen  diffusion  coefficients  (DKni)  of 
species  t,  which  depends  on  the  distance  (spacing)  between  two 
CNTs  (dcNTs).  are  of  the  same  order  of  magnitude  as  of  ordinary 
diffusion  coefficients  (Dy)  at  the  base  case  conditions. 

Similarly,  the  effect  of  Knudsen  diffusion  on  the  performance 
of  an  ordered  nanostructured  CCL  at  the  CNT  spacing  of  25  nm 
is  shown  in  Fig.  12.  It  can  be  observed  that  there  is  an  apprecia¬ 
ble  difference  in  the  performance  of  an  ordered  nanostructured 
CCL  between  the  Knudsen  and  non-Knudsen  diffusion  when  the 
CNT  spacing  is  reduced  to  half  (dcNTs  =  25  nm)  of  the  base  case  CNT 
spacing.  This  is  due  to  the  Knudsen  diffusion  coefficients  {DKn  i)  of 
species  i,  which  are  an  order  of  magnitude  smaller  than  the  ordinary 
diffusion  coefficients  (Dy )  when  the  CNT  spacing  is  reduced  to  half 
of  the  base  case  CNT  spacing  (dCNTs  =  25  nm).  Moreover,  because 
of  straight  pores  and  high  porosity  in  an  ordered  nanostructured 
CCL,  Knudsen  diffusion  effect  can  be  assumed  negligible  above  the 
CNT  spacing  (dcNTs )  of  25  nm.  Therefore,  the  effect  of  Knudsen  diffu- 


Fig.  12.  Effect  of  Knudsen  diffusion  on  the  performance  of  an  ordered  nanostruc¬ 
tured  CCL  at  dcNTs  =  25  nm. 
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Fig.  13.  Effect  of  Nation  thickness  (teflon)  on  the  performance  of  an  ordered  nanos- 
tructured  CCL. 

sion  is  not  considered  for  the  base  case  CNT  spacing  (dcNTs  =  50  nm), 
implying  the  simulation  results  presented  in  the  previous  and  sub¬ 
sequent  sections  consider  only  the  ordinary  diffusion  in  the  pore 
phase  of  an  ordered  nanostructured  CCL. 

8.  Effect  of  Nafion  thickness  ( teflon) 

Fig.  13  shows  the  effect  of  Nafion  film  thickness  (tNafion)  on  the 
performance  of  an  ordered  nanostructured  CCL.  All  the  operating 
and  design  parameters  are  kept  same  as  the  base  case  parame¬ 
ters  listed  in  Table  4.  However,  it  is  worthwhile  to  mention  that 
changing  a  design  parameter  changes  the  overall  volume  of  the 
computational  domain,  which  in  turn  change  the  volume  fractions 
of  pore  phase  and  Nafion  phase  of  an  ordered  nanostructured  CCL. 
It  can  be  seen  that  the  performance  of  an  ordered  nanostructured 
CCL  exhibits  mixed  trends  with  increasing  the  thickness  of  Nafion 
film.  The  performance  of  an  ordered  nanostructured  CCL  decreases 
with  the  increase  of  Nafion  film  thickness  for  current  densities 
up  to  1 A  cm-2.  However,  beyond  the  current  density  of  1 A  cm-2, 
the  ordered  nanostructured  CCL  with  thinnest  Nafion  film  shows 
poorer  performance.  This  is  due  to  trade-off  between  the  dissolved 
oxygen  transport  and  the  proton  transport  with  increasing  Nafion 
film  thickness,  and  the  variations  of  activation  overpotential  and 
ohmic  overpotential  with  current  density  at  different  Nafion  film 


Fig.  15.  Distributions  of  dissolved  oxygen  concentration  at  the  reaction  bound¬ 
ary  along  the  length  of  CNT  in  an  ordered  nanostructured  CCL  at  different  Nafion 
thicknesses  at  ij0  =  0.7  V. 

thicknesses  in  an  ordered  nanostructured  CCL  is  shown  in  Fig.  14. 
It  can  be  observed  from  Fig.  14  that  increasing  the  Nafion  film 
thickness  in  an  ordered  nanostructured  CCL  increases  the  cath¬ 
ode  activation  overpotential;  on  the  other  hand,  the  cathode  ohmic 
overpotential  decreases  with  the  increases  of  Nafion  film  thickness. 
This  is  due  to  increase  in  resistance  to  the  diffusion  of  dissolved 
oxygen  to  the  reaction  boundary  with  increasing  Nafion  film  thick¬ 
ness,  resulting  in  the  reduction  of  dissolved  oxygen  concentration 
at  the  reaction  boundary,  as  can  be  evident  from  Fig.  15.  However, 
decreasing  the  Nafion  film  thickness  results  in  the  reduction  of 
Nafion  volume  fraction  in  an  ordered  nanostructured  CCL,  which 
in  turn  reduces  the  effective  conductivity  of  proton  in  the  Nafion 
film,  and  thereby  increases  the  ohmic  overpotential  with  decreas¬ 
ing  Nafion  film  thickness.  The  optimal  thickness  for  the  Nafion  film 
in  an  ordered  nanostructured  CCL  is  found  to  be  lOnm  from  the 
present  simulation,  which  is  in  accordance  with  the  Nafion  film 
thickness  reported  by  Middelman  [3]. 

9.  Effect  of  CNT  length  (LCnt) 

The  effect  of  different  CNT  lengths  on  the  performance  of 
an  ordered  nanostructured  CCL  is  shown  in  Fig.  16.  The  other 
operating  and  design  parameters  are  kept  similar  to  the  base  case 


Fig.  14.  Variations  of  activation  overpotential  (a)  and  ohmic  overpotential  (b)  in  an  ordered  nanostructured  CCL  with  current  density  at  different  Nafion  film  thicknesses. 
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Fig.  16.  Effect  of  CNT  length  (LCnt)  on  the  performance  of  an  ordered  nanostructured 
CCL. 


parameters  listed  in  Table  4.  Again,  the  ordered  nanostructured  CCL 
performance  exhibits  mixed  trends  with  increasing  CNT  lengths. 
It  can  be  observed  that  increasing  the  CNT  length  increases  the 
performance  of  an  ordered  nanostructured  CCL  below  the  current 
density  of  1.5  A  cm-2.  However,  beyond  the  current  density  of 
2  A  cm-2,  the  ordered  nanostructured  CCL  with  3  |xm  CNT  length 
exhibits  better  performance  than  the  base  case  CNT  length  (5  pan). 
This  is  again  due  to  trade-off  between  the  activation  overpotential 
and  the  ohmic  overpotential  of  an  ordered  nanostructured  CCL 
with  increasing  CNT  thickness,  and  the  variations  of  activation  and 
ohmic  overpotentials  with  current  density  at  different  CNT  lengths 
is  shown  in  Fig.  17.  It  can  be  seen  from  Fig.  17  that  increasing  CNT 
length  in  an  ordered  nanostructured  CCL  decreases  the  activation 
overpotential  and  increases  the  ohmic  overpotential.  This  can 
be  attributed  to  more  reaction  sites  with  increasing  CNT  length, 
resulting  in  the  reduction  of  activation  overpotential;  in  contrary, 
increasing  the  CNT  length  increases  the  transport  distance  for  pro¬ 
tons,  resulting  in  increasing  ohmic  overpotential  with  increasing 
CNT  length.  The  optimal  CNT  length  for  an  ordered  nanostructured 
CCL  is  found  to  be  3  |jim  from  the  range  of  CNT  lengths  considered 
in  the  present  simulation. 


Fig.  18.  Effect  of  CNT  spacing  (dcNi)  on  the  performance  of  an  ordered  nanostruc- 
tured  CCL 


10.  Effect  of  CNT  spacing  (dcNTs) 

The  effect  of  various  CNT  spacings  on  the  performance  of  an 
ordered  nanostructured  CCL  is  shown  in  Fig.  1 8.  The  other  operating 
and  design  parameters  are  kept  similar  to  the  base  case  parame¬ 
ters  listed  in  Table  4.  Moreover,  Knudsen  diffusion  is  not  considered 
here  even  with  the  CNT  spacing  of  25  nm.  It  can  be  seen  that  increas¬ 
ing  the  CNT  spacing  results  in  a  considerable  drop  in  performance 
of  an  ordered  nanostructured  CCL.  Increasing  the  CNT  spacing  in  an 
ordered  nanostructured  CCL  decreases  the  number  of  CNTs  being 
accommodated  in  a  fixed  volume.  Consequently,  the  volume  frac¬ 
tion  of  Nation  decreases  with  increasing  CNT  spacing,  resulting  in 
the  increase  of  resistance  to  proton  transport.  As  a  result,  cathode 
ohmic  overpotential  increases  significantly  with  increasing  CNT 
spacing,  which  can  be  evident  from  Fig.  19.  Moreover,  it  can  be 
seen  from  Fig.  1 9  that  the  CNT  spacing  has  a  negligible  effect  on  the 
cathode  activation  overpotential.  This  is  due  to  insignificant  effect 
of  porosity  in  an  ordered  nanostructured  CCL,  which  increases  with 
the  increase  of  CNT  spacing.  Further,  the  CNT  spacing  of  50  nm  is 
found  to  be  optimum  for  the  range  of  CNT  spacings  considered  in 
the  present  simulation. 


Fig.  17.  Variations  of  activation  overpotential  (a)  and  ohmic  overpotential  (b)  in  an  ordered  nanostructured  CCL  with  current  density  at  different  CNT  lengths. 
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Fig.  19.  Variations  of  activation  overpotential  (a)  and  ohmic  overpotential  (b)  in  an  ordered  nanostructured  CCL  with  current  density  at  different  CNT  spacings. 


11.  Conclusions 

A  3D  mathematical  model  of  an  ordered  nanostructured  CCL  has 
been  developed  for  PEM  fuel  cells.  The  model  includes  all  the  rel¬ 
evant  processes  in  different  phases  of  an  ordered  nanostructured 
CCL.  Also,  the  model  takes  into  account  the  effect  of  Knudsen  dif¬ 
fusion.  The  model  predicts  not  only  the  performance  of  an  ordered 
nanostructured  CCL  at  various  operating  and  design  conditions  but 
also  predicts  the  distributions  of  various  fields  in  different  phases 
of  an  ordered  nanostructured  CCL.  The  predicted  nanostructured 
CCL  performance  with  estimated  membrane  overpotential  is  val¬ 
idated  with  measured  data  found  in  the  literature,  and  obtained 
a  good  agreement  between  the  model  prediction  and  measured 
result.  Additionally,  a  parametric  study  is  also  performed  to  exam¬ 
ine  the  effect  of  key  design  parameters  on  the  performance  of  an 
ordered  nanostructured  CCL.  In  the  absence  of  liquid  water,  it  is 
found  that  gaseous  oxygen  diffusion  in  the  pore  phase  is  not  the 
limiting  factor  for  the  performance  of  nanostructured  CCL,  owing 
to  its  parallel  gas  pores  and  high  porosity.  However,  diffusion  of 
dissolved  oxygen  through  the  Nation  film  has  a  significant  effect 
on  the  performance  of  nanostructured  CCL.  Further,  it  is  found  that 
the  effect  Knudsen  diffusion  is  negligible  on  the  performance  of 
an  ordered  nanostructured  CCL  above  the  CNT  spacing  (dcNTs)  of 
25  nm.  Furthermore,  the  optimal  thickness  for  the  Nation  film  in  an 
ordered  nanostructured  CCL  is  found  to  be  lOnm,  consistent  with 
the  Nation  film  thickness  reported  by  Middelman  [3]. 
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